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Abstract 



An equation to describe nearly ID traveling-waves patterns is put forward in the 
form of a dispersive generalization of the Newell- Whitehead-Segel equation. Trans- 
verse stability of plane waves is shown to be drastically altered by the dispersion. 
A necessary transverse Benjamin-Feir stability condition is obtained. If it is met, 
a quarter of the plane-wave existence band is unstable, while three quarters are 
transversely stable. Next, linear defects in the form of grain boundaries (GB's) are 
studied. An effective Burgers equation is derived, in the framework of which a GB 
is tantamount to a shock wave. Asymmetric GB's are moving at a constant velocity. 

PACS numbers ('American"): 47.27.Te; 47.52. +j; 03.40.Kf 
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1. Introduction 



The Newell- Whitehead-Segel (NSW) equation is a well-known universal equation to govern 
evolution of nearly one- dimensional (ID) nonlinear patterns produced by a finite- wavelength 
instability in isotropic two-dimensional media, a classical example being the Rayleigh-Benard 
convection0. This equation was applied to a number of problems, e.g., description of point 
defectsi and of linear ones in the form of the so-called grain boundaries (GB's)i. 

Thus far, the NWS equation was not extended to the case of traveling waves, i.e., to 
patterns produced by an oscillatory instability, a well-known example being the traveling- 
wave convection in a binary fluids (see, e.g.,@) (though, another extension was developedii, 
aimed at casting the the usual NWS equation into a rotationally covariant form). An 
objective of this work is to put forward a dispersive NWS equation for the description of 
nearly ID traveling- wave patterns, and apply it to relatively simple problems, viz., transverse 
stability of plane waves, and GB's. It will be demonstrated that the dispersive terms strongly 
alter the results in comparison with the usual NWS equation. In particular, all the plane 
waves turn out to be transversely unstable if a special condition of the Benjamin- Feir type 
is not satisfied; in the opposite case, the transverse perturbations destabilize exactly a 
quarter of the plane waves' existence band. Recall that, for the usual NWS equation, a 
well-known result is that a half of the band is transversely unstable against the so-called 
zigzag perturbations^. 

It is relevant to note that experiments with the traveling-wave convection were usually 
conducted in effectively ID geometries (narrow channels)!!. Recently, detailed experimental 
results were obtained for 2D patterns in a large-aspect-ratio traveling-wave convection cellEI 
(a theoretical analysis of the global patterns containing GB's was also presented i ill on the 
basis of the free-energy functional for the real 2D Swift-Hohenberq equation). These results 
demonstrate stable existence of, generally speaking, moving GB's separating domains with 
different orientations of the traveling waves. This experimental observation was an incentive 
to develop the analysis presented below. 
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The NWS equation can only apply to description of the special case when the waves 
have a small "refraction angle" at the GB. In the general case, the description of a quasi- 
stationary GB (the one that may be moving, but at a constant velocity) is based upon 
coupled Ginzburg-Landau (GL) equations^. Although the case comprised by the NWS 
equation is more restricted, a more far-reaching analysis will be developed for this case. 
An effective Burgers equation governing dynamics of GB's between the traveling waves will 
be derived, which allows one to describe not only a steadily moving GB, but also certain 
transient processes and interactions between parallel GB's. 

2. The Dispersive Newell- Whitehead-Segel Equation 

As a starting point, one can take a general two-dimensional equation for a complex order 
parameter u in an isotropic medium (cf., e.g., the derivation of the rotationally invariant 
extension of the usual NWS equation ini): 

u t = eu - (V 2 + if u + ifiy 2 )u - d\u\ 2 u - c 2 |Vu| 2 u - c 3 (Vm)V. (1) 

Here, e is a small overcriticality parameter, < e < 1, the operator (V 2 + l) 2 assumes, 
as usual, that the instability onset takes place (at e = 0) isotropically at the wave vectors 
k 2 = 1, the term f(V 2 )u with an arbitrary real-valued operator function /(V 2 ) accounts 
for the linear dispersion in a general form, and, under the standard assumption of the 
weak nonlinearity, the cubic terms with arbitrary complex coefficients ci, C2, and C3 account 
for nonlinear saturation of the instability and for nonlinear dispersion. The latter three 
terms exhaust all essentially independent local cubic terms in an isotropic system near the 
instability threshold!. 

Following the derivation of the classical NWS equation^, one considers solutions of the 
form u(x, y,t) = U (x, y, t) exp(ix +if(—l)t), where the envelope U is ssumed to be a slowly 
varying function in comparison with the rapidly oscillating carrier wave exp(ix + if(-l)t). 
Then, a straightforward asymptotic expansion leads to the following equation for the enve- 
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lope function: 



U t = eU- (2id x + d 2 y ) 2 U + i [f (-1) (2id x + flj) + (f (-1) - 2f(-l)) flg 



-(ci + c 2 -c 3 )|f/| 2 t/. 



(2) 



The convective term (the one with the first derivative) can be eliminated by a change of the 
reference frame, x — > x— 2/'(— Next, the small overcriticality e can be partially scaled out 
by means of the transformation x — > x/ \/e, y ^ y/ \fe. Here, one notices a drastic difference 
from the scaling transformation inherent in derivation of the classical NSW equation: the 
presence of the linear dispersion dictates to choose the same power of e in rescaling the 
longitudinal and transverse coordinates x and y, while, in the usual case, the scaling factor 
for y would be e" 1 / 4 instead of e -1 / 2 . 

Additionally rescaling the variables in an obvious way and properly redefining the coef- 
ficients, one arrives at the following final form of the dispersive NWS equation: 



U - (id x + V~edl) 2 U - (1 + ia)\U\ 2 U + i(3U xx + i-yU^, 



(3) 



where a, (3, and 7 are, respectively, real effective coefficients of the nonlinear, longitudinal 
linear, and transverse linear dispersions. It is important that, generically, these coefficients 
are ~ 1, i.e., they do not depend upon e, on the contrary to the coefficient in front of d 2 
in the linear operator in the second term on the r.h.s. of Eq. (JJ). Obviously, the disper- 
sive NWS equation, on the contrary to the usual one (including the rotationally invariant 
generalizations!!'!) does not admit a gradient representation. 

It should be stressed that Eq. @ is but a simplest model for the propagation of nearly ID 
waves in a nonlinear medium. Indeed, this equation takes into account only unidirectional 
waves, ignoring the waves which may travel in the opposite direction. One may expect that 
the waves traveling in both directions will be described by two coupled NWS equations, 
similarly to the coupled GL equations^. The corresponding analysis will be, of course, 
more complicated. Another simplification is the tacit omission of a possible additional real 
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equation for a slowly relaxing mode, which represents concentration field in the binary- 
fluid convection0 and occurs in other contexts^. It is known that interaction with this 
additional real mode can essentially alter properties of the solutions, especially their stability. 
Consideration of a generalized system, coupling the NWS equation to the slow mode, is 
defered to another work. 



3. Stability of Plane Waves 

As it is implied by the derivation of the NWS equation, an unperturbed plane wave can be 
chosen with the zero ^/-component of its wave vector: 

U (x, t) = a e lkx ~ lw \ ajj = 1 - k 2 , u = aa 2 + (3 k 2 . (4) 

The plane waves (f|) exist inside the band k 2 < 1. It is of obvious interest to find which 
solutions inside this band are stable against small perturbations. 

If one confines the stability analysis to pure longitudinal (y-independent) perturbations, 
Eq. (H) is nothing else but the familiar complex GL equation, in which stability criteria for 
the plane waves are well knowiS Therefore, it makes sense to concentrate here on con- 
sideration of pure transverse (x-independent) perturbations. This amounts to substituting 
into Eq. (D 

U(x,y,t) = e ikx ^ k2t V(y,t), (5) 
which yields the equation 

U t = U - (k- V^d 2 fu- (1 + ia)\U\ 2 U + i-iU yy . (6) 

Of course, the consideration of the weakly overcritical case, e C 1, implies k 2 1. 

Notice that, at k < 0, Eq. (^j) is the so-called complex Swift-Hohenberg equation, which 
was first introduced ini, and which has recently attracted a lot of attention as, e.g., a model 
dynamical equation for laser cavities0. 
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The standard approach to the stability problem is based on presenting a general solution 
to Eq. (|6|) in the form U = a exp(i(j)), with the separated amplitude and phase, in terms 
of which the plane-wave solution (^) takes the form a = ao, <p — —uit. Next, one replaces 
Eq. (|6]) by a system of coupled real equations for these variables and linearizes them with 
respect to small perturbations 8a and 5(f) of the amplitude and phase. Finally, one is seeking 
for perturbation eigenmodes is the formally complex form, 

5a = a ie CTt+ira , 5<j) = <j> x e'* Hm , (7) 

where a\ and are infinitesimal amplitudes of the perturbation, q is an arbitrary real wave 
number (independent of k in Eq. (^)), and a(p) is the instability growth rate sought for. 

After a straightforward algebra, one can cast the resolvability condition for the linearized 
perturbation equations into the form 

a 2 + 2a (al + 2k^~eq 2 + eq 4 

+ q 2 [e 2 q e + Ake 3/2 q 4 + (2a 2 e + Ak 2 e + 7 2 ) q 2 + Ak^eal + 2a^a 2 \ = 0. (8) 

Because Eq. (g) is a quadratic equation, the stability condition, Recr(g) < 0, is tantamount 
to demanding that the coefficient in front of the linear term and the free term must be 
non-negative at all real q. In the absence of the dispersion, the positiveness of the free 
term immediately leads to the well-known condition k > 0, whilst the solutions with k < 
are unstable against the so-called zigzag perturbations^. In the present case, however, the 
dispersion drastically alters this result. Indeed, taking into regard that e is a small parameter, 
while the coefficients a and 7 are not small, one obtains the following stability condition 
providing for the positive definiteness of the free term in Eq. (^|): 

07 > 0, (9) 

which IS, clS db matter of fact, a transverse version of the Benjamin- Fair (BF) stability con- 
dition well known for the longitudinal perturbations in the GL equationEi. 
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Proceeding to the condition for the positive definiteness of the coefficient in front of the 
linear term in Eq. (^), one sees that this is automatically satisfied if k > 0. That is why this 
condition did not produce any additional restriction in the case of the usual NWS equation. 
In the present case, both positve and negative k may be stable if the inequality holds 
(and both are unstable if it does not hold). A simple analysis shows that, at negative k, the 
coefficient in question remains positive at all values of q if k 2 < ag, or, finally, with regard 
to Eq. (1), if 



k > -1/V2. (10) 

It is easy to check that the dispersion equation (||) does not produce any additional 
stability condition but @ and fllPp. Thus, the eventual result of the stability analysis 
for the transverse perturbations is: all the plane waves are unstable if the transverse BF 
condition (|9|) does not hold; in the opposite case, the transversely unstable subband inside 
the existence band, — 1 < k < +1, is (according to Eq. ( |T0D — 1 < k < —l/y/2. Thus, while 
in the classical NWS equation the transverse perturbations render a half of the existence 
band, k < 0, unstable^, in the dispersive case they destabilize exactly its quarter (in terms 
of k 2 ), provided that the BF condition (Op is satisfied. 



4. The Grain Boundaries 

The consideration of GB patterns in two-dimensional nonlinear dissipative media has at- 
tracted a lot of attention starting from the worksS An analytical description of GB's 
between domains ("grains") occupied by plane waves with different orientations is usu- 
ally possible when the problem can be reduced to an effectively ID one!'!. In the case 
of Eq. (|3|), there are two natural ways to arrive at a ID problem: to substitute either 
U(x,y,t) = e tqy V(x,t) or Eq. (5), with an arbitrary constant wave number q or k. The 
former substitution, obviously, brings Eq. (0) into the form of the traditional complex GL 
equation. In terms of this equation, the GB will correspond to a shock-wave solution, which 
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is produced by a collision of traveling waves with different wave numbersEl'li. This solu- 
tion can be obtained in an approximate analytical form by means of the phase-diffusion 
equation^. Coming back to the interpretation of this solution as a GB in 2D patterns, 
one will have the GB nearly parallel to the "crests" of the colliding waves (recall that the 
transition from Eq. (1) to Eq. (|3|) implies separation of the slowly varying envelope and 
rapidly varying carrier wave whose wave vector is parallel to the x axis). The linear defect 
of this type should be, actually, described as an array of dislocations! rather than a smooth 
GB. 

The substitution (|5]), which transforms Eq. (H) into Eq. leads to essentially new 
results for smooth GB's. In this case, one should expect a GB exactly or nearly perpendicular 
to the wave crests, with a small "refraction angle" at the GB, see Fig. 5 ini. To describe 
the GB of this type in an analytical form, it is natural to consider, first of all, the lowest- 
order approximation, omitting the small terms (which vanish along with e) in Eq. (^]). This 
approximation leads to a simple but, nevertheless, nontrivial equation, viz., a special form 
of the complex GL equation without the diffusion term: 

U t = (1- k 2 )U - (1 + ia)\U\ 2 U + ijUy y . (11) 



In terms of Eq. (0), the GB again corresponds to a shock wave. As well as in the 
previous section, one should substitute U = a exp(i</>), to transform Eq. ([□]) into the 
coupled real equations, 

a t — (1 — k 2 )a — a 3 — 2^a y q — jaq y , (12) 
a(pt = —era 3 + ja yy — -yaq 2 , (13) 

where q = <p y . Following the usual assumption on which the phase-diffusion approximation 
(called the geometric optics approximation in0) is based, i.e., that the local amplitude a 
and the local wave number q are varying at a large spatial scale, one may neglect, in the 
zeroth-order approximation, the spatial derivatives on the right-hand sides of Eqs. (|T^) and 



(Pf). In this approximation, Eq. ( |12|) yields a = a/1 — k 2 . A straightforward consideration 



shows that the next-order term in Eq. (|12"D is —^aq y . Taking it into account as a small 
perturbation, one obtains the corrected expression for the amplitude: 

a= y/T=k2-l(l-k 2 )~ 1/2 q y . (14) 
This should be inserted into Eq. (|T3|) . Prior to this, it is convenient to transform Eq. 



(13), omitting the last term which is of a higher order of smallness, dividing it by a, and 
differentiating the resultant equation in y, so that to eliminate the phase by means of the 
identity <p ty = q t . This procedure produces the equation 

q t = —2aaa y — 2jqq y . (15) 

At last, substitution of Eq. ( JTS ) into Eq. ( |T3 ) yields, in the lowest nontrivial approximation, 
exactly the classical Burgers equation^ 

Qt = ( a l)Qyy ~ 2 lQQy (16) 

Notice that this equation is well-posed if the above stability condition (Q) holds. It is 
noteworthy that the combination (1 — k 2 ) does not show up in Eq. (|T6|). 

Before proceeding to consideration of its relevant solutions, it is also interesting to con- 
sider an extended equation produced by adding higher-order terms to Eq. (|i~6|). First, one 
can restore the small terms ~ y/e and ~ e, that were omitted when deriving Eq. flTT| ) from 
Eq. (^j. Because these terms are already small, it will be sufficient to calculate them at 
the lowest order in the sense of the phase-diffusion equation. Next, the second term on the 
right-hand side of Eq. (0) produces the next-order phase-diffusion correction to Eq. (|14|) . 
Collecting all the corrections, Eq. (14) becomes 

a = VT^- k 2 Y 1/2 q y - (l - k 2 Y 1/2 V~eq 2 {k + ^q 2 ) 

-\(l-k 2 y V \q y y (17) 

Finally, insertion of this equation into Eq. (|T5j ) leads to the following extended Burgers 
equation: 

q t = (<yy)q yy - 2 (7 - 2a^fe k^j qq y 



+ 4eaq 3 q y - 7 2 (l - k 2 ) 1 a {qq yy ) y + ~q y q yy + q yyy 



(18) 



The Burgers equation (|T6| ) can be exactly solved by means of the Cole-Hopf 
transformations. In particular, its simplest nontrivial solution describes a steadily mov- 
ing shock wave: 

q(.V, t) = qi + 7- 7-377 — TTT' ( 19 ) 

1 + exp (a 1 (q 2 - qi)(y - ct)) 

where q\ and q 2 are the asymptotic values of the local wave number at y = ±00 (for the 
definiteness, it is assumed q 2 > qi), and the shock wave's velocity is 

c = 7(gi+<? 2 ). (20) 

To interpret this solution in terms of the underlying 2D traveling-waves pattern, one 
should recollect the expression for the full local wave vector implied in the derivation of Eqs. 
(|) and (|TT|). According to Eq. (|5|), it is (1 + k,q), with sufficiently small k and q. This 
means that the "crests" of the waves which are forming the pattern are nearly parallel to the 
y axis, while the GB itself is strictly perpendicular to this axis. The corresponding picture 
is essentially tantamount to Fig. 5 ini. The asymptotic angles between the crests and the 
axis y (the normal to the GB) at y — > ±00 are, obviously, ~ q^ 2 ■ If Qi — —£[2, the pattern 
is symmetric, the GB being its symmetry axis, and, according to Eq. (|20|) , its velocity is 
zero. An asymmetric GB is possible as well, moving at a nonzero velocity given by Eq. (p0|). 
The same general inference was formulated ini, where GB's between traveling waves with 
nonsmall "refraction angles" were considered within the framework of coupled complex GL 
equations. It is also interesting to compare this situation with that described by the usual 
NWS equatiton. For that case, it was demonstrated ini (section III.E) that a stationary GB 
was possible only in the strictly symmetric case (qi = q 2 ), and only under a special condition 
that the asymptotic wave numbers lay exactly at the boundary of the zigzag instability, i.e., 
at k = —qf 2 - As ft was demonstrated in this section, a drastic difference of the GB between 
the traveling waves is that it exists, within the framework of applicability of Eq. (^), as a 
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generic solution. This qualitative result is also in agreement with what was obtained ini 
for the GB's with nonsmall refraction angles. It is relevant to mention that, in models of 
nonoscillatory media, a stationary GB with a nonsmall refraction angle may be asymmetric, 
but the corresponding asymptotic wave numbers must be exactly equal to onei. 

Coming back to the effective Burgers equation (fR]), one can notice that the integrability 
of this equation allows to obtain not only the shock- wave solution flUf ), but also explicit 
solutions describing transients and interactions between the shock waves, or in terms of Eq. 

interactions of the parallel GB's. In particular, solutions describing formation of a GB 
out of various initial perturbations, as well as a collision of two shock waves and their merger 
into a single one, are available0. Although a more detailed comparison with experimental 
observations of the traveling waves in the large-aspect-ratio binary-fluid convection cells@ 
is necessary, it appears that formation of a GB and merger of GB's as a result of their 
interaction can be observed. 

It may be interesting to analyze influence of the small additional terms in the extended 
Burgers equation flTB| ) on properties of the shock waves, but this issue is left beyond the 
framework of this work. Lastly, it is relevant to mention that the special GL equation ( |TTD 
has an altogether different physical application in terms of nonlinear fiber opticsEl). If t 
and y are interpreted, respectively, as the propagation distance and the so-called reduced 
time in a fiber, this is a propagation equation for the light in the presence of a broadband 
gain and linear and nonlinear lossesEi. The condition (|9[) is then equivalent to the condition 
that the fiber's dispersion is normal. In this case, cw (continuous-wave) states with different 
frequencies q are stable, and the shock wave solutions, e.g., ([19]), describing collisions between 
the cw's, are of a certain physical interest. 
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